% This is a script for reading in the phantom files, as an example
Nx = 284; Ny = 411; Nz = 722;

% Load in optical and acoustic phantoms
fid = fopen('./MergedPhantom.DAT', 'r');
phan = fread(fid, 'uint8=>uint8'); phan = reshape(phan, [Nx, Ny, Nz]);
fclose(fid);



